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1 Introduction 

Sub-Riemannian geometry deals with tiie study of smootli manifolds M that are en- 
dowed with a vector distribution A and a smoothly varying positive definite quadratic 
form. The distribution zi is a subbundle of tangent bundle TM and the quadratic form 
allows measuring distance between any two points pi,P2 G M [I], [2], [3]. Other names 
that appear in literature for sub-Riemannian Geometry are Carnot-Caratheodory ge- 
ometry [4], Non-holonomic Riemannian geometry [5] and Singular Riemannian geome- 
try [6]. The aim of defining and solving a sub-Riemannian problem is to find the opti- 
mal curves between two given points pi, p2 on the sub-Riemannian manifold M such 
that sub-Riemannian distance between the points is minimized [2 ,[3\. Sub-Riemannian 
problems occur widely in nature [2], [7] and have therefore been extensively studied via 
geometric control methods on various Lie groups such as the Heisenberg group [8], 
S^, SL(2), SU(2) ^, SE(2) [10,, Engel group [11 , Solvable groups il2j, SOLV' [S], 
and also in |14|.|15).|16|.|17|. Few examples of physical systems that describe sub- 
Riemannian problems and on which Geometric control methods have been successfully 
applied include parking of cars, rolling bodies on a plane without sliding, motion plan- 
ning and control of robots, satellites, vision, quantum mechanical systems and even 
finance [1|,[7]. 

We consider the sub-Riemannian problem on the group of motions of Pseudo- 
Euclidean plane which is a subspace of pseudo Euclidean space. The pseudo Euclidean 
space F^"^ is (n -|- m)-dimensional space defined over field of real numbers R and 
endowed with a non-degenerate indefinite quadratic form q |18] : 



An important example of pseudo Euclidean Space is the Minkowski space-time that 
arises in the special theory of relativity. It is essentially a four dimensional pseudo 
Euclidean space F^^^, with three ordinary dimensions of space and another inter- 
mingled dimension of time |19) . Minkowski space-time was a refomulation of special 
theory of relativity and it presented a mathematical setting in which Einstein's theory 
of relativity and Lorentz geometry could be mathematically formulated. 

A pseudo Euclidean plane is a two dimensional subspace F^~^^ with q{x) = x\ — x\ 
|18| . As complex numbers represent vectors on a Cartesian/Euclidean plane, hypercom- 
plex or split complex numbers are used to represent gyrovectors on pseudo Euclidean 
plane [20] . A pseudo Euclidean plane represents Mikowskian space-time of two dimen- 
sions with one spatial variable and one temporal variable [20j and is therefore associated 
to our understanding of the world. Due to its natural linkage with the Minkowskian 
space-time plane, formulation of a sub-Riemannian problem on pseudo Euclidean plane 
can possibly give insight into how nature works. 

The motions of pseudo Euclidean plane described in Section [2] form a 3-dimensional 
Lie group known as special hyperbolic group SH(2) [21]. The optimal control problem 
comprises a system of left invariant vector fields with 2 dimensional linear control 
input and energy cost functional. The group SH(2) gives one of the Thurston's three- 
dimensional geometries [22' and the study of Sub-Riemannian problem on SH(2) bears 
significance in the program of complete study of all left-invariant Sub Riemannian 
problems on three dimensional Lie groups following the classification in terms of the 
basic differential invariants [23] . 
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Notice that an equivalent sub-Riemannian problem was considered in |13| on Lie 
group SOLV~ . However, the parametrization of sub-Riemannian geodesies obtained 
in |13| is far from complete. This paper seeks rigorous scheme of analysis developed in 
|10|.|11).|16| for parametrization and qualitative analysis of the extremal trajectories. 
The corresponding results are insightful, in simpler form owing primarily to the use of 
simpler elliptic coordinates and allow further analysis on global and local optimality 
of geodesies. The paper is organized as follows. We begin with a brief description of 
SH(2) in Section 2 and introduce the concepts related to sub-Riemannian geometry in 
Section 3. Section 4 and 5 contain the main results of this research. In Section 4 we 
state the sub-Riemannian problem on SH(2) and investigate the global controllability 
of the control distribution. We apply the Pontryagin Maximum Principle on SH(2) 
and discuss the abnormal and normal trajectories. Section 5 is a detailed description 
of the integration of the vertical and horizontal subsystem in Elliptic coordinates. In 
Section 6 we present the qualitative analysis of projections of extremal trajectories on 
xy-pl&ne. Section 7 and 8 pertain to future work and conclusion respectively. 

2 The Group SH(2) of Motions of Pseudo Euclidean Plane 

Following exposition is motivated from 21 and is presented here for the sake of com- 
pleteness. 

2.1 Pseudo Euclidean Plane 

A pseudo Euclidean plane is a 2-dimensional real linear space endowed with an indefi- 
nite bilinear form given as: 

[a, b] = aibi — a2h2- 

The distance r between a point a(ai,a2) and another point b(6i,62) on the plane is 
given by the formula: 

= (ai — hif' — (a2 — 62)^ = [a — b, a — b]. 

Unlike the Euclidean plane, the distance in pseudo Euclidean plane may be real as 
well as purely imaginary. The distance between two distinct points can also be zero 
provided a\ — hi = ±(02 — &2)- Geometrically, pseudo Euclidean plane is represented 
by unit hyperbola a\ — a\ = \ that also depicts Minkowski space-time |20| . see Figure 
[1] The asymptotic lines ai = 02 and ai = —02 segregate the plane into four distinct 
sectors known as Right (RS), Left (LS), Up (US) and Down (DS) sectors. We will 
consider our problem only on RS where — ai < a2 < ai and distance between points 
r > because the motions of pseudo Euclidean plane we consider are sector preserving 
maps. As the points in a Euclidean plane are transformed into polar coordinates via 
trignometric functions, the points in RS may be represented in polar coordinates by 
hyperbolic functions as follows[20].|21|: 



ai = rcoship, 
a2 = rsmhtp, 
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Fig. 1 Pseudo Euclidean plane represented by unit hyperbola 

where r G is the length and </3 G R is the hyperbohc angle of rotation of the 
gyrovector. 

2.2 Group SH(2) of Motions of Pseudo Euchdean Plane 

The motions we consider are non homogeneous, linear, distance, orientation and sector 
preserving maps of points in Pseudo Euclidean plane. Specifically the motions comprise 
translations and hyperbolic rotations given as: 

h\ = a\ cosh z + a2 sinh z + 
62 = 11 sinh 2 + 02 cosh z -\- y, 

where x,y,z £ R. 

Thus motion m of pseudo Euclidean plane is completely parametrized by x,y,z £ 
R. Composition of two motions mi{xi,yi, zi) and m2{x2,y2, Z2) is another motion 
m3{x3,y3,Z3) given as: 

'm3{x3,y3,Z3) = rni{xi,yi,zi).m2{x2,y2,Z2), 

where, 

X3 = X2 cosh zi + y2 sinh zi + xi, 
y3 = X2 sinh zi + y2 cosh zi+yi, 

Z3 = 21 + Z2- 

The identity motion mj^ is given hy x = y = z = 0, and inverse of a motion m(x, y, z) 
is given by {x^ , z^) where. 
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= —X cosh z + y s'mh z, 
= X sinh z — y cosh z, 



1 

z 



The composition of motion m with m ^ is given as m = mj^. 

2.3 Lie Group and Lie Algebra Representation 
The group SH(2) can be represented by third order matrices: 

{/ cosh z sinh 2 a; \ "] 
sinh z cosh z y \ x,y,z & R\ . 
\ 1] J 

The Lie group SH(2) comprises three basis one-parameter subgroups given as: 

(cosht sinht 0\ / 10 t\ 

sinhf cosht , W2{t) =010, W3{t) = 
1/ \0 1/ 

whereas basis for Lie algebra are the tangent matrices = \t=o to the sub- 

groups of Lie group SH(2). Ai are given as: 



Ai=100 , yl2 = 000L A 







The Lie algebra is thus: 

£ = I/rfM = sh(2) = span{Ai,A2, A3} . 

The multiplication rule for C is [^4, B] = AB — BA. Therefore, the Lie bracket for sh(2) 
is given as [^1,^2] = A3, [^1,^3] = A2 and [A2, A3] = 0. 



3 Sub-Riemannian Geometry 

3.1 Riemannian Manifold 

Smooth Riemannian space/manifold {M,g) is a real smooth manifold M on which 
an inner product Qp can be defined on each tangent space TpM, Vp G AI. The inner 
product varies smoothly everywhere on AI such that for any vector fields X and Y 
on AI and p £ AI, p 1— > gp{X{p),Y{p)) is a smooth function [2],[^. The family of all 
inner products gp on AI allows to define various geometric constructs such as length 
of curves, distance between points and angles on Riemannian manifold exactly like the 
scalar product defines such notions on Euclidean space. This family of inner products 
is termed as Riemannian metric (tensor) [24j . 
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3.2 Sub-Riemannian Manifold 

Sub-Riemannian space/manifold is a variation of the Riemannian manifold. It com- 
prises a manifold M of dimension n, a smooth vector distribution A whose rank m is 
constant such that m < n, and g is a Riemannian metric on A. It is denoted as a triple 
(M, A, g). The distribution Z\ on M is a smooth linear subbundle of the tangent bundle 
TM i.e. A C TM. Intuitively, on sub-Riemannian manifold, the motion is restricted 
along paths that are tangent to horizontal subspaces or the admissible directions of 
motion are constrained to horizontal subspaces Aq, q £ M |25|.|26|. Sub-Riemannian 
manifolds naturally arise in such diverse areas as non-holonomic systems in classical 
mechanics, image reconstruction, image inpainting etc see e.g. [2] .[7] .1271.1281. 



3.3 Sub-Riemannian Distance 

Consider a sub-Riemannian manifold (M, A, g) and a Lipschitzian horizontal curve 
7 : J C R — >■ M; ^(t) G ^-y(t) for almost all t £ I. The length of 7 is given as: 



length{'y) = J ^ g^(t){i{i))di, 
I 

where g-y(t) is the inner product in zi..y(()0. The sub-Riemannian distance between 
two points p,q £ M is length of the shortest curve joining p to q: 

, ■ r \ , , / \ 1 is horizontal curve \ 
d(p;g)=zn/|/enfftM7): ^ Jo^r^s p to q 



3.4 Sub-Riemannian Problem 

Consider a driftless dynamical system on sub-Riemannian manifold {M,A,g): 

m 

q = ^Ui{t)fi{q), gM™. 

i=l 

The problem of finding horizontal curves 7 from initial state go to final state qi with 
shortest sub-Riemannian distance ci(go; ?i) and tangent to a given distribution Aq C 
TqM is called a sub-Riemannian problem [^.[271. Intuitively there exists a set of vector 
fields fi whose values ^q G M form a local orthonormal frame of the sub-Riemannian 
structure {A; g). The horizontal curves 7 : / C R — > Af are the solutions of the following 
optimal control problem in M: 

m 

q = J2'^di)Xt{q)), q€M, (ui , ■ ■ ■ , G R™, 

i=l 

q{0)=ld, q{ti) = qi, 



J 
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4 Sub-Riemannian Problem on SH(2) 

Consider following driftless control system on SH(2): 

q = Uifi{q)+U2f2{q), qeM = SH(2), (ui.ua) G R^ (1) 
q{0) = ld, q{ti) = qi, (2) 

I = j + u\ dt — ^ min, (3) 

fj_{q) = qA2, f2{q)=qAi. (4) 

Here ([TJ represents the dynamical system with control inputs Ui and control distri- 
bution A = span{/i,/2}. In ([2|, g(0) is the intital state at time t = and 
represents the final state to be reached at time ti whereas / is the sub-Riemannian 
distance (length functional) to be minimized. Canonical frame on AI in terms of |23| 
is given as: 

/2(g), fo{q)=qA3, 

[/l,/0]=0, [/2,/o]=/l, [/2,/i]=/o. (5) 

By |23j . the sub-Riemannian structure: 

iM,A,g), A = span{fi,f2}, g{.h, fj) = 5ij , 

is unique upto rescaling, left invariant contact sub-Riemannian structure on SH(2). 
Here Sij is the Kronecker delta. 

Hypothetically, the system represents a Reeds-Shepp car moving on a hyperbolic 
plane. The car is a rigid body and is allowed to move forward, backward and rotate 
about a vertical axis passing through its centre of mass. Being point mass, the position 
and orientation of car on hyperbolic plane are the coordinates (x, y) of its centre of mass 
and the angle z respectively. Given an initial state qo and a final state qi , the objective 
is to determine a horizontal curve 7 C SH(2) between the initial and final state that 
minimizes the sub-Riemannian length functional (|3|. The admissible trajectories are 
assumed Lipschitz and admissible control inputs are assumed bounded. 



4.1 Complete Controllability 

Theorem 1 The control system ^ is completely controllable. Moreover, the optimal 
control problem ^^-^ has solutions. 

Proof Consider the family of vector fields A C Vec{Sii{2)), given by ([l]), where 
yec(SH(2)) is the set of all vector fields on SH(2). Since [fi{q) , /2(g)]= -[.f2(q), fi{q)] 
= —fo{q), the system is full rank because the distribution A satisfies the bracket gen- 
erating condition (also known as Hormander condition). Hence, 

£,Z\ = sparz{/i(g),/2(g),-/o(g)} =TqSH(2) Vg £ M. 

By Rachevsky-Chow's Theorem |29|.|30|. for a connected manifold and corresponding 
bracket generating control distribution, the system is completely controllable. From 
the definition of SH(2) it is connected and A is bracket generating, hence the system 
([1]) is completely controllable. 

Existence of optimal trajectories for problem ([l])-(|4| follows from Filippov's theo- 
rem [25]. □ 
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4.2 Pontryagin's Maximum Principle for Sub-Riemannian Problem on SH(2) 
In coordinates (i, y, z) the basis vector field are given as: 

fi{q) = coshz^ + sinhz^, 

and 

Therefore ^ may be written as: 







/ cosh z \ 






■ 


H 


sinh z 


|.t+| 








{ 







(6) 



By Cauchy Schwarz inequality, 

/ ti \ ^ ti 



{l{u)f = / ^uj+u^dt <ti {ul + ul)dt 



Thus sub-Riemannian length functional minimization problem (|3| is equivalent to the 
problem of minimizing the following energy functional with fixed ti ^26j: 

*i 

1 f 2 2 

J = - / (til + U2)dt — >• min. (7) 


We write the PMP form for ([l},([2|,([7| using coordinate free approach described in |25| . 
Consider control dependent Hamiltonian for PMP corresponding to vector fields /i (q) 
and /2(g): 

hU^) = {\,fu{q))+'^{ui+ul), q = n{\), XgT*M. (8) 

Let hi(X) = (A, fi{q)) be the Hamiltonians corresponding to basis vector fields /j. Then 
Q can be written as: 

K{X) = uihi{X) + U2h2{\)+'^{ul+u^), ueR^. (9) 

Now PMP for optimal control problem is given by using Theorem 12.3 |25) as: 

Theorem 2 Let u{t) be optimal control and q(t) be optimal trajectory for t € [0, ti] 
and hu{X) given by ^ be the Hamiltonian function for Then, there exists 

a nontrivial pair: 

iu,Xt)j^O, u€K, XteT*^(t)M, 7v{Xt) = q{t). 

where At is a Lipschitzian curve and v G { — 1,0} is a number, for which following 
conditions hold for almost all time t £ [0,ti]; 

At = ^fi(t)(At), (10) 
^fi(t)(At) = max/i;;(4)(Af). (11) 

where h'^f^^.^^Xt) is the Hamiltonian vector field corresponding to the Hamiltonian func- 
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4.3 Abnormal Trajectories 

Abnormal trajectories correspond to the case u = 0. The Hamiltonian ^ in this case 
can be written as: 

h°{\)=uihi{X) + U2h2{\). (12) 

Theorem 3 All abnormal extremal trajectories for problem ([T| , ([2| , ([7| are constant. 

Proof Let At be the abnormal extremal, then the maximization condition of PMP 
yields: 

hi{\t) = h2{\t) = 0. (13) 



Differentiating (131 w.r.t. Hamiltonian vector field and noting that Poisson bracket 
follows the same multiplication rule as that of Lie bracket of sh(2): 

hi = /iij = {uihi + U2h2,hi} = ui {hi, hi} + M2 {/i2, ^i} = U2ho, 

h2 = h2^ = {uihi + U2h2,h2} = ui {hi,h2} + «2 {h2, ^2} = — mi/iq- 

Therefore, 

U2{t)ho{\t) = Ul{t)ho{\t) = 0, 

ul{t)ho + wi(f)/io = 0. 

If /lo(At) = for some At, then /ii(At) = /i2(At) = ho{Xt) = which means At = 0. 
This is impossible since 1/ = 0. Therefore Ui{t) + u^it) = iti = U2 = and 

hence the abnormal extremal trajectories are constant. □ 



4.4 Normal Trajectories 

Normal trajectories correspond to the case v = —1. The Hamiltonian Q in this case 
can be written as: 



H 



hZ\\)=Ulhl{\)+U2h2{\)~'^{ul+ut^ , ueV?. (14) 



Using the maximization condition of PMP, the trajectories of the normal Hamiltonian 
satisfy the equalities: 

dH _ / fti — iti 

du \h2 — U2 



0, 



ui=hi, U2 = h2- (15) 

Note that if Ui = 0, then normal extremal trajectories are constant. Therefore the 
abnormal trajectories are not strictly abnormal. The normal extremals are the trajec- 
tories of Hamiltonian system A = -fz (A), A £ T*M with the maximized Hamiltonian 
H = \ {hl + hi) > 0. Specifically, for non-constant normal extremals H > 0. Note 
that Hamiltonian function in normal case is homogeneous w.r.t. hi, /12 and therefore 
we consider its trajectories for the level surface H = ^. The phase cylinder containing 
the initial covector A in this case is: 



C = TlMn{H{X) = ^1 = {{hi,h2,ho) €R^\hl + hl = l}. (16) 
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Differentiating (141 w.r.t. Hamiltonian vector field we get: 



hi = {H,hi} = 1^ (^/if + hf^ = h2{h2,hi} = 

/i2 = {H, 112} = I ^ (j^i + ' '^aj = hi {/ii, /12} = -hiho, 
ho = {H, /lo} = I ^ (^hl + hfj ,hQ^ = hi {hi, ho} + ^2 {^2,^0} = hih2 
Hence, complete Hamiltonian system in normal case is given as: 



fhl\ 


/ h2ho \ 


h2 


-hiho 


ho 


hih2 


X 


hi cosli z 


y 


hi sinli z 


\ i / 


\ h2 ) 



(17) 



Theorem 4 Vertical subsystem of the Hamiltonian system in normal case is a 
mathematical pendulum. 



Proof Introduce following coordinates transformation: 

hi = cos a, /i2 = sin a. 

Thus, 

hi = — sin aa = sin q./iq, 
a = ~ho. 



Similarly, 



ho = cos a sin a = ^ sin 2a 



Let us introduce another change of coordinates: 

7 = 2q G 25'^ = R/AtvZ, c = -2ho G M, 

=^ 7 = 2ti = ~2ho = c, 



and 



Thus, 



c = ~2ho = — 2/ii/i2 = — 2 cos a sin a = — sin 2q = — sin 7. 



c 

- sin 7 



(18) 



(19) 



(20) 



(21) 



It can be easily seen that ( 22 1 represents a mathematical pendulum. 



(22) 
□ 
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5 Parametrization of Extremal Trajectories 



5.1 Hamiltonian System 



Hamiltonian system for normal trajectories was given in (17 1. Under the transforma- 
tions introduced in (1181), 1211), the horizontal subsystem can be written as: 







/ hi cosh z \ 




■ 




hi sinh 2 













cos ^ cosh z ^ 
cos ^ sinh z 
sin I , 



(23) 



5.2 Decomposition of the Initial Phase Cylinder 

Following the techniques employed in |10| . the decomposition of phase cylinder C pro- 
ceeds as follows. The total energy integral of the pendulum obtained in ( 22 I is given 
as: 

2 

E = — -cos'y = 2hl- hl + hi, £;g[-1,+oo). (24) 



The initial phase cylinder ( 16 1 may be decomposed into following subsets based upon 
the pendulum energy that correspond to various pendulum trajectories: 



5 

c=[jc\, 

where 



Ci 


= {A G C\E e 


(-1,1)}, 


C2 


= {A G C\E G 


(1,(X))}, 


C3 


= {A G C\E = 


l,c^O}, 


Ci 


= {A G C\E = 


-1} = {(7,c) G Cl7 = 27rn,c = 0}}, n G N, 


Cs 


= {A G C\E = 


1} = {(7,c) G C|7 = 27rn-|-7r,c = 0}, nGN 



Continuing the approach taken in [lOj the subsets Ci may be further decomposed as: 



Ci = ul^oCl, C\ = {(7,c) G Cils<7n(cos(7/2)) = (-1)*}, 
C2 = C+ U C^, Ct = {(7,c) G C2\sgnc = ±1}, 

C3 =uLo(C^+UCr), ={(7,c)GC3jsffn(cos(7/2)) = (-l)\sffnc = ±l}, 

C4 = ul^oCl, Cl = {(7, c) G Cl7 = c = 0}, 
Cs = yj]=oCl, Cl = {(7, c) G C|7 = 27Ti + 7T,c = 0}. 



In all of the above i = 0, 1. The phase portrait of the pendulum and corresponding 
decomposition of initial phase cylinder C is depicted in Figure |2] 
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5.3 Elliptic Coordinates 



Employing the approach developed in |10).|11) we transform the system in terms of el- 
liptic coordinates {(p, k) on the domain U^_^Ci C C. Note that ip is the reparametrized 
time of motion and k is the reparametrized energy of the pendulum. Correspondingly, 
we describe the system and the extremal trajectories in terms of Jacobi elliptic func- 
tions sn(^, fc), cn{ip,k), dn{ip,k), a.m{ip,k), and E{(fi,k) = J^^ dn'^(t,k)dt. Detailed 
description of Jacobi elliptic functions may be found in |31) . 



5.3.1 Case 1 : \ = {ip,k) G Ci 



■ 7 

sm 2 = •'I 



1 
2 
c 



ksn{(p,k), si = sgn ^cos ^ j , 
dn(^, fc), 
kcn{(p,k), <p £[0,4K{k)]. 



cos ^ = sidn(^, k) 



(25) 
(26) 
(27) 
(28) 



Proposition 1 In elliptic coordinates the flow of vertical subsystem rectifies. 



Proof Using ( 25 1 



,2 • 2 7 , c' 
k = sm ^ + 



.2 

2 ' 4 ■ 



Taking the time derivative of ( 29 1 



Using (221 



2fcfc = 2sin|cos^^ + | 



2kk = ^ sin 7 — ^ sin 7 = 0. 



(29) 



(30) 
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Now either fc = or fc = 0. Since k £ (0, 1), hence it cannot be zero and therefore: 

k = 0. (31) 



Using (281 and the derivatives of elhptic functions defined in |31| . 
d 



d rc 
dt V2 



dt 



kcn(ip, k), 



c , d , .-.dip , d , , dk , , ^ dk 



2 dip dt dk 

■ sin7 = — 2fcsn(iy3, fc)dn(<^, k)ip 



dt ' 



because ^ = 0. Now, 



Now using (261,1271 



sin7 

2fc sn((/p, fc) .dn((/p, k) 



sin — cos — = siksn{ip, k).sidn{ip, k), 



(32) 



Thus ( 32 1 becomes: 



T T 2 
2 sin — cos — = 2siksn{ip, k).dn{ip, k), 

sin7 = 2k sn{ip, k).dn{ip, k). 
if = \. 



(33) 
□ 



5.5.a Case 2 : \ = {ip,k) € C2 



sm 



E+1 



S2sn I pfc 



cos — = cn I — 



1 



■ 2 7 , 

sm^ J + X 
S2 = sgn(c). 



G (0,1), 



'^^-dnl^.fc), ip€[0,'lkK{k)]. 



(34) 

(35) 
(36) 
(37) 



5.3.3 Case 3 : X = {ip,k) G C3 



k = l, 

sin ^ = S1S2 tanh si = sgn ^cos ^ j , S2 = sgn(c), 

7 

cos ^ = COSh(/3, 

- = S2/ COSh(/3, ip £ ( — CXD, 00). 



(38) 
(39) 
(40) 
(41) 



Using the procedure outlined for Case 1, it can be proved that the flow of the pendulum 
rectifies for cases 2 and 3 as wefl. 
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5.4 Integration of Vertical Subsystem 

Since the fiow of vertical subsystem rectifies in elliptic coordinates, therefore, the ver- 
tical subsystem is trivially integrated as ip = t + (fio and k = constant, where ipo is the 
value of (/5 at f = 0. 



5.5 Integration of Horizontal Subsystem 



In the following we consider integration of horizontal subsystem ( 23 1 for cases 1-3 noted 
above. Assuming zero initial state i.e. a;(0) = y{0) = z(0) = 0. 



5.5.1 Case 1 : \ = {(p,k) G Ci 

Theorem 5 In case 1 extremal trajectories are parametrized as fo 



illows: 



2 



V 



si In [(dn (fi — ken ip) .w] 



where w = -i -, . 



(42) 



written as: 



Proof From (23 1 consider z = sin ^ = sifc sn((/?, fc). The solution to this ODE can be 

(43) 



z H> 
J dz = J siksmpdif. 
fo 



Using [32j, eq(43l becomes 



z = si ln(dn ip — ken ip) — si ln(dn ^ ken tpo)- 
Let \nw = -ln(dn(po -fecni^^o), w = yo~fccn yo ' Then (44 1 becomes 

z = si ln[(dn ip — fccn ip).w]. 



From (231 now consider. 



X = cos — cosh z = 



si 

X = — 



sidn<^cosh (si In [(dn;^ — kcn(p).w]) , 

■^i ^ 2 , , , duifi \ 

— w.dn tp — kw.dn pen p + j- n — 

2 Y (dn<^ — Kcn<^).TOy 



(44) 



(45) 



(46) 



This can be integrated 



as: 



si 

X = — 



V 



f , 2 , , /" 1 7 1 /" dn 05 -|- fccn <z)dn p 

w / dn pap — kw / dn (pen ipdip H / — - — = t-f. — 7^ 

I J w J dn-'ip — k-'en-'d 



fo 



'ip 



dp 



(47) 



Title Suppressed Due to Excessive Length 
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Now using the standard identities of elliptic functions result of integration of ( 47 1 can 
be written as: 



si 



w + 



1 



.(l-fc2) 



[E(^) - E(^o)] + 



k 



w(l~ fc2) 



— kw ] [sn </p — sn (po] 



(48) 



From (|23|) now consider 
, 7 



y = cos — sinh z = sidn (p sinh(si ln[(dn (p — ken (p).w]), 

y = sfdn ip sinh(ln[(dn p — ken p).w]), 
y = dn ip sinh(ln[(dn ip — ken 'p).w]). 



(49) 



The integration follows the same pattern as described above and hence final result of 



integration of ( 49 1 can be written as: 
1 



y ■ 



w - [E{p) - E(^o)] - +kw][snp- sn^o] 



(50) 
□ 



5.5.2 Case 2 : \ = {<p,k) G C2 



Theorem 6 Consider horizontal system for case 2 (34) -(31) and substitute 
ipo = "jf and = J = 4>o + ^, the integration results can be summarized as: 

1 / 1 



2 \w{\ - fc2) 
1 /, k 



S2 



1 



2 \lt)(l-fc2) 
S2 I , k 



where w ■ 



z = S2 ln[(dn t/j — fccn iIj).w\ 
1 



E(^)-E(^o)-fc'(V'-V'o) 
[snV" - sni/)o] , 
+ w ] [E(V) - E(^/>o) - k'^{ilj - Vo) 
antf) — sntpo] , 



dn ipQ—kcn ^0 

Proof Proof follows from the procedure outlined in case 1. 



(51) 



□ 



5.5.3 Case 3 : \ = {<fi,k) G C3 

Theorem 7 In this case, k = 1. Integration results are summarized as: 



j^{ip- Pq) + w (tanh ip - tanh po) 
^ {(p ~ ifio) - w (tanh ip - tanh ipo) 
— S1S2 \n[w sechp] 



(52) 



and w = cosh ipo . 
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Proof Consider horizontal system (23 1 for case 3 (38l-(41l 



i = sin — = S1S2 tanh (p, 

z = — S1S2 [ln(sech (/p) — ln(sech</3o)]- 

Let — ln(sech </Po) = In it), it) = coshiyso, then: 

2 = — S1S2 ln[it) sech 99]. 

From ([23| now consider, 

i = cos — cosh 2 = sisech ip cosh (— S1S2 ln[i;; sech , 

■ SlSech(^ r ln[twscchip] I — In [tw scch 

X - - + e 



dx 



sisech ip 



1 + w'^sech^ip 
w sech <p 



dip, 



2 



1 



((/p — (/pq ) + TO (tanh p — tanh po ) 



(53) 



(54) 



From (231 now consider. 



y = cos — sinh 2 = sisech p> sinh (— S1S2 ln[it) sech p]) . 



y 

dy 



S2 sech p ^^ln[w sech ^] ^— ln['u; sech ip] j 

-S2 sech 



S2 

2 



; sech p — [w sech lys] ^\ dp 



{p — po) — ui(tanh p — tanh po) 



(55) 
□ 



5.6 Integration of Horizontal Subsystem - Degenerate Cases 

In the following we present the integration of horizontal subsystem in degenerate cases 
i.e. A G C4 and A G C5. 



5.6.1 Case 4 : X € C4 

Theorem 8 Integration results in case 4 are summarized as follows: 

sgn (cos 2 ) ^ \ 


/ 



(56) 



Title Suppressed Due to Excessive Length 



17 



Proof 

Since z{Q) = 0, 
Therefore, 



^ = sin-=sm(^— j =0. 



z = Q. (57) 



a: = cos ^ cosh z = cos ^^^^ > 



i = sgn {cos | j , 

X = sgn {cos ^) * + Wx, 

X = sgn {cos t, (58) 
where Wx = becasue a;(0) = 0. Now, 

y = cos ^ sinhz = cos {^~^~^ sinh(O) = 0, 

y = Wy 

y = Q, (59) 
where Wy = because y{0) =0. □ 

5.6.2 Case 5 : \ € C5 

Theorem 9 Integration results is case 5 are summarized as follows: 

' \ 

. (60) 

sgn (sin J 




Proof 



Thus, 



^. 7 _ . / TT + 2n7r 



2; = sm - = sm I ) = sgn 



(sin 2) 



z = sgn (sin^j t + Wz, 

z = sgn {sin | j t, (61) 
where Wz—O because z{0) = 0. Now, 

X = cos % cosh z = cos I ^ "^^"^ j (,Qgij ;j: = 0, 



2 

a: = 0, (62) 

because a;(0) = 0. Now, 

y = cos ^ sinh z = cos { ^ ] ginh z = 0, 



2 V 2 , 

y = 0, (63) 
because j/(0) = 0. □ 
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6 Qualitative Analysis of Projections of Extremal Trajectoris on xy-Flane 

The standard formula for the curvature of a plane curve (x{t),y{t)) is given as |33) : 



xy — xy 

(±2+^2)^ 



(64) 



Using (17 1,(64 1 curvature of projections {x(t),y{t)) of extremal trajectories of the 



Hamiltonian system (17 1 is given as 



2 cos^ ^ (cosh 2z) i 



(65) 



The curves have inflection points when sin 1 



and cusps when cos Z 



or c = 0. 

We see that all curves [x{t),y{t)) have inflection points for A £ uf_j^Ci but only for 
A G C2 the curves have cusps. The resulting trajectories (x{t),y{t)) are shown in 
Figures [3] |4] [5] In degenerate case 4 i.e. A G C4, the extremal trajectories qt are sub- 
Riemannian geodesies in the plane {z = 0}. The curve {x{t), y{t)) is a straight line on 
the X-axis. In case 5 i.e. A G C5, the curve {x{t), y{t)) is just the initial point (0, 0) for 
{x = y = Q}. For non-zero initial conditions x{Q) = Wx, 2/(0) = Wy, the motions of 
pseudo Euclidean plane are only hyperbolic rotations whereas the translations along 
a;-axis and j/-axis are zero. The resulting trajectory is a quarter circle in the RS that 
approaches the upper and lower arms of the hyperbola in RS of Figure [l] as t — >■ 00. 




7 Future Work 



Most natural extension of this work is the computation of Maxwell strata and obtaining 
the global bound on cut time based on discrete symmetries of the vertical subsystem. 
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Fig. 4 Inflectional Trajectories A £ C2 

y 




Fig. 5 Inflectional Trajectories A £ C3 

To this end the methods developed in |10).|17|.|16|.|14) have been employed and re- 
sults shall be reported in another paper. Investigation of local and global optimality 
of extremal trajectories via description of conjugate and cut loci is another exciting 
and challenging dimension in the problem under consideration. Description of global 
structure of exponential map and optimal synthesis is the ultimate goal to be addressed 
in the entire research on SH(2). 

8 Conclusion 

Sub-Riemannian problem on SH(2) is important from perspective of Mathematics 
as well as Physics. The direct connection between the pseudo Euclidean plane and 
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Minkowski space-time geometry suggests that such analysis could potentially lead to 
better understanding of the Special Theory of Relativity. The Mathematics perspective 
is also equally significant. The group SH(2) as abstract algebraic structure has its 
own significance and sub-Riemannian problem on SH(2) is important in the entire 
program of study of three dimensional Lie groups. In this paper we have obtained the 
complete parametrization of extremal trajectories in terms of Jacobi elliptic functions 
and described the nature of projections of extremal trajectories on xy-plane. This shall 
pave the way for further work on the goals highlighted in Section [7] and is target of our 
future work on SH(2). 
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